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ric particles at the LHC, using the new matrix element generator SUSY— MADGRAPH . The cross 
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amplitudes to fermions, and derive sum rules for weak supersymmetric couplings. 
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I. INTRODUCTION 



The matter content of the highly successful Standard Model of particle physics is gener- 
ally considered to be fully revealed after the discovery of the top quark in 1994, although 
the exact mechanism of clectroweak symmetry breaking remains undetermined [1]. The 
Standard Model description of spontaneous symmetry breaking is minimal, involving only 
one additional complex scalar doublet. This introduces an unsatisfactory instability in the 
scalar sector of the theory. A theoretically more attractive scenario is that spacetime re- 
spects the maximal extension of the Poincare symmetry, supersymmetry (SUSY) [2]. Its 
minimal version, the MSSM, simultaneously provides solutions to several problems in high 
energy physics and cosmology: a candidate for weakly-interacting dark matter; possible 
unification of the gauge couplings at high energies; and stability of the scalar sector which 
generates electroweak symmetry breaking through renormalization group running. 

SUSY must be a broken symmetry at low energy, as we do not see spin partners of the 
Standard Model particles. As a result, the squarks, sleptons, charginos, neutralinos and 
gluino of the MSSM must be massive in comparison to their Standard Model counterparts. 
Experiments such as LEP and Tevatron [3] have put stringent bounds on some of the SUSY 
partners' masses. It will fall to the LHC to perform a conclusive SUSY search covering 
masses all the way to the TeV scale. Real physics will, however, begin only after a potential 
SUSY discovery: in particular the strongly-interacting squarks and gluinos can be produced 
in large numbers at the LHC [4, 5], and their decay cascades typically carry kinematical 
information about a large fraction of the weakly interacting SUSY spectrum [6, 7]. This 
information can be used to narrow down different SUSY breaking mechanisms [8, 9]. 

Much MSSM and non-minimal SUSY phenomenology has been performed over the years 
in preparation for LHC, nearly all of it using relatively simple dominant 2^2 processes at 
leading order or next-to-leading order [4, 5]. Often, these calculations involve a number of 
approximations, many of which might not be sufficient for practical applications once the 
collection of data begins [10] . Examples include: consideration of spin correlations and finite 
width effects in SUSY particle production and decays; SUSY-electroweak (EW) and Yukawa 
interferences to some SUSY-QCD processes; exact rather than common squark masses in 
the t-channel; and additional 2 — > 3 or 2 — > 4 particle production processes, such as the 
production of additional hard jets in squark and gluino production [11], or weak boson fusion 
(WBF) production of colorless SUSY particles, the subject which we address here. 

The last topic is particularly interesting, because it may help us to observe sleptons and 
weak inos at LHC. These particles can be extremely difficult to observe in direct production 
channels due to small rates and very large backgrounds, and their appearance in squark and 
gluino cascade decays can depend on the SUSY breaking scenario. WBF is an electroweak 
process which naturally leads to high-p^ but very far forward taggable jets and little central 
jet activity. It can be used to observe small electroweak signal rates in a region of phase 
space not very populated by QCD events. This was applied very successfully to heavy [12] 
and intermediate-mass Higgs boson production in the Standard Model [13], as well as weak 
boson scattering [14], the MSSM Higgs sector [15], and to multi-doublet Higgs models [16]. 
As a consequence, there has been recent interest in using the same technique for colorless 
intermediate-mass SUSY pair production, with investigations of a few channels in limited 
regions of parameter space [17-19]. The reported results were often negative, but appeared 
to hint at some promising regions of parameter space [20]. Searching for sleptons and weak 
inos in these channels can provide useful information about the SUSY Lagrangian. 
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We present a comprehensive calculation of supersymmetric colorless WBF production 
channels at the LHC, including a broad scan over parameter space. We also discuss some 
theoretical issues regarding consistency of electroweak gauge couplings and unitarity of weak 
boson scattering to weak inos, and their relevance for practical calculations at LHC energies. 

To compute these production rates we introduce a new, supersymmetric version of the ma- 
trix element generator madgraph ll [21], which properly takes into account various physics 
aspects which are usually approximated in the literature, such as those listed above. It also 
provides a practical means to include hard jet radiation effects [11] and eventually match to 
parton shower simulations, as well as to evaluate more complicated production processes in 
LHC and linear collider phenomenology [10]. In Sec. II we discuss SUSY— MADGRAPH and 
its model assumptions, 2 — 2 unitarity as a test of its MSSM couplings, requirements for 
MSSM electroweak couplings to maintain unitarity of weak scattering processes at high en- 
ergy, and gauge invariance of WBF processes. We give results for LHC WBF pair production 
of colorless SUSY particle in Sec. HI. We present conclusions in Sec. IV and a technical 
overview of MADGRAPH n and SUSY-madgraph in the Appendix. 

II. SUSY-MADGRAPH, UNITARITY AND MSSM COUPLINGS 

The new matrix element generator susy-madgraph is an extension of the Standard 
Model MADGRAPH H package with the new feature of Majorana fermions. ^ The implemen- 
tation of MADGRAPH H into the web-based event generator madevent [21] is straightfor- 
ward and has already been used for the calculations published in Ref. [11]. To include the 
complete MSSM particle spectrum we use a pair of model data files describing the particle 
content and its interactions. Specifically, we describe the MSSM as the minimal supersym- 
metric model which conserves -R-parity, does not contain any CP- violating complex phases, 
and is CKM- and MNS-diagonal. However, those assumptions can be dropped by straight- 
forward changes in the code. We do not assume any particular SUSY breaking scheme, so the 
MSSM spectrum and couplings can be handed to it by any spectrum generator, regardless 
of what assumptions go into constructing the spectrum. 

A. Setup and Tests 

The most work-intensive step in extending MADGRAPH H to inchidc the MSSM is the 
correct definition of all couplings in terms of the Lagrangian parameters. To implement the 
MSSM couplings, we work with the conventions of Refs. [22] and Ref. [23], and cross-check 
with those of Ref. [24]. Due to the multiple couphngs conventions in the literature [22, 24- 
26], we employ a large number of numerical checks to ensure their correctness. In addition 
to a number of numerical comparisons to published SUSY cross sections at e+e^ [27] and pp 
colliders [5], we also test the MSSM couplings to electroweak bosons by checking unitarity 
of 2 — > 2 scattering processes at very high energy. This not only serves to ensure that 
the couplings are correct, but also reveals a general complication in the use of spectrum 
generator input from the SUSY Les Houches Accord (SLHA) [28]. SLHA is a standard- 



Codo available for download at littp://pheno. physics. wise. edu/"^lehn/siiiadgraph or http://www. 
pas . rochester . edu/~rain/ smadgraph 
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ized format for communicating SUSY Lagrangian and low-energy parameters, from spec- 
trum generators [29, 30] and SUSY particle widths [31] to multi-purpose Monte Carlos and 
next-to-leading order predictions [5] . susy-madgraph uses this convention for the SUSY 
spectrum input. We discuss this complication in subsection II C. Finally, we check some 
500 2 — > 2 cross sections numerically [10] between SUSY-madgraph and the multi-purpose 
event generators WHIZARD [32] and SHERPA [33] . 



B. General VV —>■ xx Unitarity Sum Rules 



A powerful analytical and numerical check of couplings which enter the production of two 
fermions in gauge boson scattering are unitarity sum rules, which describe the scattering 
process in the hmit of large center-of-mass energy E. The general process reads: 



(1) 



The incoming and outgoing particle masses are mi^2 and Mi 2, respectively. The incoming 
gauge boson polarizations are Ai_2, and the final-state fermion helicities are ai^2 = L, R. 
Four types of Feynman diagrams can contribute to the above process: ^-channel exchange of 
a fermion of mass M^. ?i-channcl exchange of a fermion Fi with mass M^, annihilation to 
an s-channel vector boson V, or to an s-channel scalar S. The helicity amplitudes for these 
four diagrams comprise the matrix element TW^IJ^^) froni which we derive the sum rules. 

The amplitude AiQ(^ can be written in terms of general couplings g±. For exam- 
ple, the interactions between two fermions Fi, F^ and a gauge boson V appear in the 
Lagrangian as C = g^^^'^^ Fij^PaF2V^. The left- and right-handed couplings a = ± 
correspond to the left- and right-handed projectors P-t = (1 ± 75)/2. Similarly, fermion- 
scalar couplings appear as £ = g±^'^^ FiP±F2S, triple-vector-boson couplings as £ = 

and vector-vector-scalar couphngs as £ = gViV2S3 ViV2nS. 

We take z to be the incoming beam direction and define the outgoing particles to lie in 
the x — z plane. The scattering angle sin^ describes the fraction of the final state momentum 
in the x direction: 
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The remainder represents terms which do not increase with E. Each line in Eq. (2) 
corresponds to one class of diagrams. The vector boson Vs may be 7 or Z. We note that 
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the amplitude Eq.(2) and all equations in this section are obtained in the unitary gauge. A 
scalar field S, therefore, represents both the Higgs bosons and the Goldstone bosons. To 
obtain sum rules in the Feynman gauge, remove all terms proportional to — m\. 
Prom Eq. 2 we can derive two sum rules for J = and J = 1: 
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Similarly, the amphtude M^q is given by 
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leading to two more sum rules 

J = 0: Yl [2M,g'-'''''9? 



F2V2 



Migl'^'^'gl""^'^' - M2gl'"'''g'J 

-,FiFiV2 ^FiFiVi _ F1F1V2 F1F2V1 



Fk 



FiFaVg* 



V3 



M^gF.F2Vz'^ m| 



(6) 



5 



J= 1 : 



J2 [-Mi^^^^^^^y^^^^^^ - M^/^^*^^/^^^^^] 

Fe 

+ E [{-Migr-''' - ) gv.v.vj ^ . (7) 



The mixed-helicity amplitudes are 
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and lead to two sum rules 
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We note, however, that the last two sum rules are not independent from the first four. 
They can be derived from the hehcity-diagonal cases in the hmits Mi^2 = 0. 

To illustrate how the sum rules work we give their explicit form for the process W~Z — > 
XiXi- The rule for {R,R) and J = in terms of general couplings is given in Eq. (3). If 
we assign Vi = W~ ,¥2 = Z and Fi — Xi , F2 — , we can exchange i-channel neutralinos 
Fk = Xk M-channel charginos Fg = xj ■ The sum over s-channel scalars does not 
appear in this example because the charged Higgs docs not couple to the initial-state gauge 
bosons. The couplings g± for this case are given in the appendix. The sum rule we have to 
numerically check now becomes: 



y-2m^ig^ +^xr^^ 3l +^x°,9r 9r 
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To verify these SUSY— madgraph couplings we numerically check the set of sum rules for 
the processes W'^W" — > xtxj: W''Z — > Xi)^ ~^ XiXs- The last process also serves 
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as a check for the proper description of neutrahnos with negative mass eigenvalue. In that 
case we can either use a phase to re-rotate the mass matrix onto positive eigenvalues, which 
means working with a complex mixing matrix, or we can use the negative mass eigenvalue 
and make use of an analytic continuation of the expression for the matrix element [27]. The 
two approaches are equivalent at leading order and respect the unitarity sum rules. 

In addition to these sum rules, we check unitarity numerically for amplitudes produced 
by SUSY-MADGRAPH. Our test includes more than 300 (2 — > 2) scattering processes with the 
initial states VV,V'~f and VHi, where V = W^, Z and — h, H, A, H^; and the final states 
Xi'^Xj'"^'^ 7 ^ Hi, ^Hi, HiHj and /i/j • Wc vary the center-of-mass energy from threshold to 
IC^ TcV, to avoid problems with machine precision. Wc require that amplitudes at most 
approach a constant at high energy. Unfortunately, 2 — > 2 unitarity is not sufficient to check 
triple-Higgs (or any triple-scalar) couplings. Instead, we verify them by comparison with 
pubhshed results [34]. 

C. Unitarity and the use of SLHA input 

Let us consider a (2 — >• 2) scattering process VV — >• X^XpXtxJ ^ ^ = W^,Z at high 
energy. There is a well-known gauge cancellation between s-, t- and u-channel diagrams, 
just as in any Standard Model process VV — > //. This limits the high-energy behavior of 
the amplitude to at most approach a constant. For the cancellation to occur, all couplings 
and masses in the scattering must be exactly related by a consistent set of electroweak 
parameters, driven by gauge invariance of the Lagrangian. Even the slightest deviation from 
this condition will result in the amplitude growing as E or E^, depending on the process. 
This provides a powerful test of the neutralino and chargino couplings. However, when we 
check the high-energy behavior for weak ino combinations using SLHA input from a default 
SUSY spectrum generator, the test fails. The answer lies in an inconsistent treatment of 
electroweak parameters between the default assumption of most spectrum generators, and 
the way collider processes are conventionally calculated. 

The neutralino and chargino sector is described by the mass matrices 

^ 


-mzSyjCp nizCyjCfs 
\ mzSy,S/3 -mzCyjSfs 

with the bino, wino and higgsino mass parameters on the diagonal and the gaugino-higgsino 
mixing masses in the off-diagonal elements. The neutralino and chargino mixing matrices 
which enter the matrix element are computed with a given set of electroweak parameters 
mw, mz and Sw{= sin 9^)- The same weak parameters also enter the matrix element through 
the Standard Model gauge boson couplings. Both sets must be consistent to assure proper 
cancellation of the different diagrams at high energy. This is similar to Standard Model 
top quark pair production, where the final state masses and the Yukawa coupling have to 
be identical. The reason for the test failure is because the electroweak parameters used for 
collider calculations are generally taken to be those at the Z pole, while most spectrum gen- 
erators run the electroweak parameters up to some scale Q to diagonalize the mass matrices, 
where the default scale Q is often the SUSY scale, 0{1) TeV. Even though the electroweak 
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FIG. 1: Cross sections with the inconsistent electroweak parameters relative to that with the 
consistent SLHA ripping scheme described in the text, as function of the hadronic center-of-mass 
energy. The red (sohd) curve is for xtxi pairs, the blue (dashed) curve for xtxi production, and 
the green (dot-dashed) curve for The cuts of Eqs. (14) and (15) were imposed on the two jets. 

parameter differences at these two scales may seem to be quite small, the mismatch violates 
gauge invariance, i.e. scattering unitarity is not preserved. 

Because the renormalization group evolution implemented in spectrum generators pre- 
dicts DR parameters at the TeV scale, our solution to preserve unitarity is simply to extract 
the electroweak parameters from the neutralino and chargino mass matrices given in the 
SLHA input. From the form of the mass matrices in Eq. (13), extracting the effective elec- 
troweak parameters is straightforward. If these matrices are loop-improved, this approach 
is equivalent to assuming a (yet to be explored) universahty [35] and absorbing the loop 
corrections into the effective electroweak parameters mf^,m|^,s^ and tanf]^^. We fix the 
overall coupling strength hj Gp = 1.16639 x 10~^ GeV~^, and extract mz, rnw and from 
the mass matrices. If the weak mixing angle is defined in the on-shell scheme, these three 
parameters will be consistent also at the TeV scale. These effective parameters are obviously 
not the pole masses or proper MS masses. Therefore, we should examine how large an effect 
running of the electroweak parameters can have on practical calculations at colliders such 
as the LHC. 

We quantify the impact of the different electroweak parameter schemes in Fig. 1, where 
we show the relative rates of three particular WBF chargino and neutralino production 
cross sections with and without the SLHA ripping scheme described above. For the SUSY 
spectrum, we select the generic SPSla parameter point. At realistic collider energies the 
technical gauge-invariance violation leads to cross section deviations at a fairly constant level 
of C(10%) for LHC through VLHC (200 TeV) energies. The fiat behavior indicates that 
we do not yet see unitarity violation numerically at LHC or VLHC energies. Nevertheless, 
(9(10%) corrections from the scheme change exhaust the typical QCD error bars to WBF 
cross sections [36]. 
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FIG. 2: Cross section with varying Yukawa coupling relative to that with cro(^/ = ^/) for WBF 
Tifi (red solid) and hib*i (blue dashed) production. The cuts of Eqs. (14) and (15) are imposed 
on the two jets. The crosses represent the points where the Yukawa coupling is extracted from the 
masses and as described in the text. 

Similar to the above-described problem with electroweak gauge invariance and unitarity, 
another issue arises in the squark-Higgs and slepton-Higgs couplings. The symmetric scalar 
mass matrix is completely defined by three entries, usually chosen as two mass eigenvalues 
and the mixing angle. They can be computed from the left- and right-handed soft-breaking 
masses, the quark Yukawa coupling, the trilinear mass parameter Aq and the Higgs-sector 
parameters /i and tan/?. However, for example in the f f*hP coupling the same parameters 
mj, y4j, yU, tan/5 appear in a different combination from the off-diagonal mass matrix entry. 
If we compute the matrix element for the production process WW —>■ titj, these Lagrangian 
parameters appear in the couplings to an s-channel Higgs, but they also enter implicitly 
through the stop mixing angle (given fixed stop masses). Again, there is potential for 
a mismatch in the matrix element. In contrast to the electroweak parameter mismatch 
described above, three-scalar couplings cannot spoil unitarity in 2 — > 2 scattering processes. 
Although the mismatch breaks SUSY between the top and stop couplings, since the SUSY 
violation in the three-scalar coupling is soft, this does not reintroduce a quadratic Higgs 
mass divergence at higher loop order. 

There is again a way around this, but it is not as straightforward as in the neu- 
tralino/chargino sector. We must choose one SUSY parameter to be extracted from the 
squark and slepton mass matrices, assuming fixed mass eigenvalues. Because the stau and 
stop mass matrices are not necessarily evaluated at the same scale, the universal parameters 
tan/? and /i should not be redefined as effective parameters to cure the mismatch in the 
stau and stop sectors simultaneously. On the other hand, the renormalization of Af has a 
strong impact on the perturbative convergence of the light Higgs mass mh,_wh.i\e attempting 
to preserve SUSY by replacing the Yukawa coupling as it enters the ff*h coupling also 
requires a shift in the Yukawa coupling f fh. Strictly speaking, these couplings are fixed 
by gauge invariance. However, we also know that, taking Higgs production and decay as 
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an example, calculations using the running MS Yukawa couplings significantly improves the 
perturbative behavior. Given these complications, and given that keeping the mismatch 
corresponds to only a small shift of a soft-breaking Lagrangian parameter, we use the usual 
DR parameters m/, Aj, /x, tan/3 in the SUSY-madgraph couplings ff{h,H,A}. On the 
other hand, for practical purposes it will probably be preferable to use the running Yukawa 
coupling in the ffh couping as well as in its ff*h counterpart. 

We investigate the practical impact of this inconsistency via the rate change of WBF 
sfermion pair production (even though the rates arc too small to be observed). For stau 
pairs, the rate change is typically at the per mille level, because of the small tau Yukawa 
coupling. As another check we construct matrix elements for WBF bibl production, even 
though this is a colored final state. The rate changes for the SPS points vary from negligible 
(SPS 2,6) to almost a factor of 4 (SPS 4). We show the cross section rates with varying 
Yukawa coupling relative to those with Yf — mf in Fig. 2 for LHC f^f{' and bibl production 
at SPS la. We see that the Yukawa couphng extracted from the masses and from ^/ is indeed 
the expected DR or MS running mass at the high scale given by the SUSY masses. 

D. Interplay of WBF and non-WBF processes 

As an example, we discuss neutralino pair production, with representative Feynman dia- 
grams shown in Fig. 3. The VFVF-fusion component proceeds via both a i-channel diagram 
involving a mixed gaugino and higgsino W'^XiXf coupling, and s-channel Z and Higgs bo- 
son diagrams mediated via a higgsino- induced Zx^x^- coupling. The heavy Higgs diagrams 
essentially do not contribute, as the HVV coupling is proportional to cos(/3 — a), which 
vanishes rapidly for large Higgs masses niA ^ 160 GeV. The t-channel diagrams of the ZZ- 
fusion processes are induced by the higgsino content alone, and have only Higgs bosons in 
the s-channel. 

There is a class of bremsstrahlung diagrams where the incoming quarks exchange an 
electrowcak gauge boson, and a quark line emits a Z boson which decays into a neutralino 
pair: Fig. 3(c). One would naively estimate this contribution to be small (after the jet cuts), 
but there are gauge cancellations between the diagrams which cannot be neglected. 

There also exist diagrams where a pair of incoming quarks exchange a gauge boson in the 
t-channcl, and one quark line splits to a neutralino and a squark, which decay to neutralino 
plus quark. Fig. 3(d). In other words, this is a double-bremsstrahlung contribution, induced 
by the gaugino content of the neutralino. However, for the SPS scenarios, the squarks are 
heavier than the weak inos, so these diagrams constitute higher-order clcctroweak resonant 
neutralino-squark production. Moreover, this contribution does not produce the typical 
forward jets with central neutralinos, and is a separate gauge set of diagrams. Therefore, 
we do not include this contribution. Strictly speaking, for production of light neutralinos in 
scenarios with relatively light squarks (beyond the SPS benchmark points) these diagrams 
can be non-negligible. 

For chargino pair production, the story is much the same, although ,s-channcl Z boson 
diagrams now also have a gaugino component. For slepton/sneutrino production the WW 
and ZZ fusion diagrams as well as the Z-bremsstrahlung are identical to the neutralino 
case, but the double-bremsstrahlung diagrams do not exist. 

As a check of the consistent generation of amplitudes by susy-madgraph , we test the 
U{1)em (electromagnetic) gauge invariance of WBF production of MSSM particles. To test 
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(a) 




(b) 




(c) 




FIG. 3: Representative Feynman diagrams for WBF neutralino pair production at a hadron col- 
lider: (a) t-channel chargino exchange; (b) s-channel Z boson (or Higgs bosons, not shown); (c) 
Z bremsstrahlung; (d) double-neutralino bremsstrahlung, which in most cases is a higher-order 
electroweak correction to on-shell squark- neutralino production. 



electromagnetic gauge invariance, we generate matrix elements for WBF production plus a 
photon for a number of weak ino and sfermion pairs. These are (2 5) processes, which 
would be impossible to calculate without automated tools such as susy-madgraph. We 
then set the photon polarization vector numerically to its four-momentum, and simultane- 
ously set all particle widths to zero, to confirm that the matrix element vanishes. 

Although we had hoped that the EM gauge invariance of the generated amplitudes with 
one additional photon emission would demonstrate the need to include non-WBF amplitudes, 
mainly because the photon is partly an electroweak (14^^) boson, we found that the WBF and 
non-WBF amplitudes give rise to separately EM gauge invariant sets in all the processes 
we studied. It is interesting that there are further U{X)em gauge-invariant subsets, some 
of which are the set of all diagrams involving a four-point WSS coupling - that is, of 
two gauge bosons and two MSSM scalars - as well as diagrams involving two scalars and 
s-channel Higgs propagators, as they arise from the SUSY Lagrangian D-terms, which are 
separately gauge invariant. 

Checking only the U{X)em gauge-invariance of the WBF processes is, therefore, not a 
verification that calculations neglecting bremsstrahlung diagrams is gauge invariant with 
respect to the full electroweak theory. For this one has to check the exact electroweak Ward 
identities for each process. This is not as simple as replacing the vector boson wave function 
with its momenta, as it is in the photon case. Unfortunately, testing the full electroweak 
gauge invariance of WBF processes is beyond the scope of this paper, as it requires further 
development of automated matrix element tools. 

While previous calculations [17-19] are therefore V{V)em gauge- invariant, some of their 
results are incorrect due to large cancellations between diagram subsets, most notably for 
slepton-sneutrino production. We quantify this in the next Section. 
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III. WEAK BOSON FUSION CROSS SECTION RESULTS 



A characteristic of WBF is that the incoming quarks are scattered at relatively small 
angles, but at typical transverse momentum of order ~ My/2. These scattered partons can 
easily be identified as jets in CMS [37] and ATLAS [38] for pr > 20 GeV. We thus impose 
generic observability kinematic cuts on the (parton-level) scattered quarks: 

PT(j)>20GeV, |r/(j)|<5.0. (14) 

Unless the WBF process involves an internal photon propagator, which occurs only for 
charged particle production, almost all the cross section will pass these cuts. A further 
characteristic of WBF is that the two jets he in opposite hemispheres, i.e. far forward and 
far backward, with large pseudorapidity separation between them. We can therefore impose 
the typical WBF separation cut 

|r/(ji)-r/(j2)| >4.2, (15) 

which corresponds to a pseudorapidity separation of 3 units between the jet cone edges, 
using a typical LHC jet cones size of radius 0.6. This cut typically suppresses the WBF 
rate by a factor 3-4. SUSY-QCD backgrounds which include a pair of hard jets can be a 
few orders of magnitude larger than the WBF cross section before cuts. However, they are 
typically suppressed by 2-3 orders of magnitude by the WBF cuts. Normally, for a WBF 
analysis one imposes further constraints to ensure the WBF-produced decay products lie 
between the two jets, and the two forward jets form a large invariant mass. These cuts 
usually result in little additional loss of the signal. Because they are not necessary for our 
illustrative purposes here, we ignore such refinement of the signal rates. Forward-tagged 
jets are typically identified with these cuts with an efficiency of etag = 0.86, predicted by 
detector simulation, for low-luminosity running. This number might decrease somewhat for 
high-luminosity running [39]. 

We use leading-order CTEQ6L1 parton densities [40] in our calculations, with factor- 
ization scale fJ'f = {fTii + 1712) /2, where mi and m2 are the masses of the heavy particle 
pair produced. The electroweak parameters except for Gp are taken from the SLHA model 
files, as described in the previous section. We investigate all SPS points [41], MSSM pa- 
rameter space benchmarks designed to represent a number of canonical scenarios to align 
experimental studies and phenomenology. We consider them to be only starting points for 
investigation, with no assumption that the scenarios are more likely to be realized in nature 
than any other. 

A. Charginos and neutralinos 

Neutralinos and charginos can be pair-produced in WBF, for their rates we always impose 
the tagging jet cuts of Eqs.(14) and (15). Most cross sections arc at the attobarn [ab] 
level, with a few exceptions. Any cross section below 1 fb is almost certainly going to be 
unobservable, unless it has a large branching ratio to a particularly distinctive final state 
that can be observed with high detection efficiency, such as multiple electrons or muons. 
Even then it may have sizable backgrounds from cascade decays of colored SUSY particles 
produced at QCD strength. We comment on possible noteworthy cases. 
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We first show the neutraUno pair production cross sections in Tabic I, which are almost 
universally at the few-attobarn level. In SPS scenarios, the lightest neutralino is mostly 
bino, which reduces the coupling to weak bosons, while the higgsinos are heavy. The single 
exception to tiny neutralino pair cross sections in all SPS scenarios is X2X2 production, 
where the X2 mostly wino with a sizable higgsino fraction. For this case, rates typically 
range from a tenth of a femtobarn to half a femtobarn. At SPS9, the X2 decays mainly to 
i~^i~Xi{^ = e, II, r), which results in a sizable detection efficiency of due to the 4-lepton final 
state with large missing transverse momentum. Since there are very few Standard Model 
backgrounds which could produce this signature, SUSY backgrounds are the dominant source 
and likely small. There is some chance this channel could be observed, although with very 
low statistics unless one considers the LHC luminosity upgrade, SLHC [39]. 

However, as we show in Table V, the Drell-Yan (DY) rate for X2X2 production is far 
larger [5]. Unfortunately, the DY and WBF production mechanisms are not very comple- 
mentary, in the sense that both processes involve the same gaugino-induced coupling to 
light-flavor quarks and the higgsino-induced coupling to gauge bosons. The similar-size rate 
for WBF XiXi production at SPS9 is likely uninteresting, as the WBF Z ^ vv rate is much 
larger [42]. 

We next show the chargino pair production cross sections in Table II. For opposite-sign 
production, where we do not apply the forward tagging jet cuts as the signature is already 
highly distinctive, they are at most around a femtobarn for xXX\ SPSla and SPS9. 
However, given that the WBF W~^W~ background is several orders of magnitude larger, 
we do not expect this signal to be observable. Moreover, WBF production of opposite-sign 
charginos is again much smaller than the direct DY production. Much more interesting are 
the same-sign chargino rates, also shown in Table II. While the signal rates stay below a 
femtobarn even for xtxtt the Standard Model WBF W^W^ background is now similar 
in size to the signal. Instead, the main background would more likely come from SUSY 
cascades, for example from gluino pairs or squark pairs which decay to like-sign dileptons. 
Those production cross sections are typically above ~ 100 pb [5]. but the branching fractions 
which lead to dileptons are usually tiny. If like-sign chargino pairs were clearly identifiable 
at LHC, it would provide conclusive evidence that the neutralinos appearing in the t-channel 
are Majorana fermions, although this would first require separate evidence that the chargino 
or the neutralino candidates are in fact fermions. We are not aware of any other way of 
demonstrating this at the LHC, and extraction of this information out of cascade decays 
would likely be tedious. Of course, a future linear collider operating in e~e~ mode would 
provide definitive proof of the Majorana nature of the neutralinos. 

A complete study of signal and backgrounds for same-sign chargino production is on 
the way. The analysis breaks down into two dominant regions of parameter space for the 
signal: where the chargino decays predominantly to W^x^ followed by a leptonic W boson 
decay; and where it decays preferentially to f^Vj. and the stau ultimately decays to a lepton 
and the lightest neutralino. In both cases, the x\ gives rise to significant missing transverse 
energy. The final state is then a pair of far-forward/backward high-p^^ jets, a pair of same-sign 
central leptons (also expected to be high-pr due to the large mass of the cascade parent), and 
significant missing transverse energy. There would also be very little additional jet activity, 
due to the clectroweak nature of the scattering process, as mentioned in the Introduction. 
The efficiency for such a signal would be extremely high. The Standard Model WBF ly^iy^ 
background is also of C(l) fb [14], but should display markedly different lepton kinematics. 
SUSY backgrounds would arise primarily from QCD squark and gluino production, which are 
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often collectively at the tens to hundreds of picobarn level. However, only a small fraction of 
these will give high-p^ forward jets from the cascade decay to jets plus charginos, and much 
of the sample will have additional high-p^ jets which can be vetoed. The lepton kinematics 
will again be significantly different, as the typical relatively light charginos (compared to 
the squarks and gluinos) will yield much higher-p^^ leptons. 

Finally, we show the mixed chargino-neutralino production cross sections in Tables III 
and IV. These are also at the attobarn level, the only exception being xfxt production, 
with a typical rate of 0{1) fb. However, this should be compared to the DY rate of many 
hundreds of femtobarn as shown in Table V. 
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x?x^ 


0.018 


0.001 
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0.003 


0.008 
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0.001 
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0.001 





0.001 


0.001 


0.001 


0.002 


x^x^ 


0.52 


0.10 


0.24 


0.10 


0.26 


0.29 


0.039 


0.057 


0.15 





x^x^ 


0.049 


0.008 


0.009 


0.008 


0.026 


0.008 


0.009 


0.017 


0.016 





x^x^ 


0.065 


0.011 


0.011 


0.011 


0.034 


0.009 


0.023 


0.045 


0.022 





x^x^ 


0.006 


0.001 


0.001 


0.001 


0.004 


0.001 


0.004 


0.009 


0.003 





xlxl 


0.008 


0.001 


0.001 


0.001 


0.004 


0.001 


0.005 


0.013 


0.003 





xlxl 


0.007 


0.001 


0.001 


0.001 


0.004 


0.001 


0.008 


0.020 


0.003 






TABLE I: Cross sections [fb] for WBF neutralino pair production at LHC, for all MSSM benchmark 
SPS points, using the kinematic cuts of Eqs. (14) and (15). Cross sections are shown to two 
significant digits or rounded to the nearest attobarn, and those smaller than half an attobarn are 
shown as zero. 
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x^xr 


1.6 


0.26 


0.63 


0.27 


0.74 


0.77 


0.13 


0.23 


0.42 


1.3 


X1X2 


0.056 


0.010 


0.011 


0.010 


0.029 


0.010 


0.015 


0.028 


0.019 


0.002 


xtx2 


0.035 


0.007 


0.004 


0.006 


0.020 


0.003 


0.030 


0.068 


0.015 





xtxt 


0.93 


0.22 


0.48 


0.23 


0.51 


0.57 


0.067 


0.077 


0.31 


0.88 


x^x^ 


0.13 


0.022 


0.028 


0.022 


0.070 


0.015 


0.072 


0.14 


0.049 


0.002 


x^x^ 


0.001 











0.001 





0.011 


0.032 


0.001 





xrxr 


0.28 


0.056 


0.13 


0.058 


0.14 


0.16 


0.017 


0.020 


0.083 


0.25 


xrx^ 


0.040 


0.006 


0.005 


0.006 


0.021 


0.005 


0.018 


0.036 


0.014 


0.001 


X2X2 




















0.002 


0.007 









TABLE IL Cross sections [fb] for WBF opposite-sign and same-sign chargino pair production 
at LHC, for all MSSM benchmark SPS points, using the kinematic cuts of Eqs. (14) and (15) 
for opposite-sign charginos only (see text). Cross sections are shown to two significant digits or 
rounded to the nearest attobarn, and those smaller than half an attobarn are shown as zero. 
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0.001 
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0.001 
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0.007 
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0.001 
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0.007 
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0.006 
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0.022 


0.015 


0.001 


X^Xi 
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0.001 





0.001 


0.001 


0.001 


0.001 


X^X2 


0.029 


0.005 


0.005 


0.005 


0.015 


0.004 


0.010 


0.021 


0.010 





X^X3 


0.007 


0.001 


0.001 


0.001 


0.004 


0.001 


0.005 


0.011 


0.003 





X^X4 


0.008 


0.001 


0.001 


0.001 


0.004 


0.001 


0.008 


0.021 


0.003 






TABLE III: Cross sections [fb] for WBF positive-sign chargino plus neutralino pair production at 
LHC, for all MSSM benchmark SPS points, using the kinematic cuts of Eqs. (14) and (15). Cross 
sections are shown to two significant digits or rounded to the nearest attobarn, and those smaller 
than half an attobarn are shown as zero. 
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0.009 


0.001 


0.001 


0.001 


0.002 


0.001 


0.002 


0.006 


0.002 


0.32 


xrx2 


0.39 
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0.15 


0.062 


0.18 


0.19 


0.026 


0.044 


0.097 





xrxs 


0.027 


0.004 


0.005 


0.004 


0.014 


0.005 


0.006 


0.012 


0.009 


0.001 


xrx4 


0.027 


0.005 


0.005 


0.004 


0.015 


0.004 


0.007 


0.014 


0.009 


0.001 


X2'Xl 


0.002 











0.001 








0.001 


0.001 


0.001 


X^X2 


0.018 


0.003 


0.003 


0.003 


0.010 


0.003 


0.006 


0.013 


0.006 





X^X3 


0.004 


0.001 





0.001 


0.002 


0.001 


0.003 


0.007 


0.002 





X2X4 


0.005 


0.001 





0.001 


0.003 





0.005 


0.012 


0.002 






TABLE IV: Cross sections [fb] for WBF negative-sign chargino plus neutralino pair production at 
LHC, for all MSSM benchmark SPS points, using the kinematic cuts of Eqs. (14) and (15). Cross 
sections are shown to two significant digits or rounded to the nearest attobarn, and those smaller 
than half an attobarn are shown as zero. 





SPS la 


SPSS 




DY 




DY 


A\'BF 


x^x^' 


25.4 


0.52 


2.90 


0.15 


Xi^Xi 


705 


1.6 


212 


0.42 


X^X2 


828 


0.64 


276 


0.17 


xrx2 


474 


1.26 


142 


0.31 



TABLE V: Drcll-Yan (DY) and WBF pair production cross sections [fb] for two representative 
MSSM benchmark SPS points. The DY results were calculated with Prospino2.0 [5] at NLO. 

B. Sleptons and sneutrinos 

The second type of supersymmetric WBF processes is the production of charged sleptons 
and sneutrinos, again with the tagging jet cuts of Eqs. (14) and (15). Smuon cross sections 
are identical to those for selectrons, as their masses are identical, so we omit them in our 
tables. 

Sneutrino pair production rates, shown in Table VI, are at the attobarn level and univer- 



15 



sally too small to be observed. The charged slepton cross sections and the mixed cases (cf. 
Tables VII and VIII) show similarly limited promise (in contradiction with Refs. [18, 19]). 
At least part of this discrepancy is due to large cancellations between the pure-WBF di- 
agrams and bremsstrahlung, where the incoming quarks scatter off each other via a weak 
boson, and the sparticle pair arises from Z/W^-Bremsstrahlung off a quark line. The cancel- 
lation is practically nil for the neutral current [ZZ fusion) in slepton pair production, as well 
as sneutrino pairs, but a little more than a factor 3 in the charged current (VFIV fusion), 
which dominates, leaving an overall cancellation of about a factor 3. The cancellation is 
much larger for mixed sneutrino-slepton pairs, about a factor 30. Refs. [18, 19] considered 
only the WBF diagrams and so did not notice the cancellations. We emphasize that these 
cancellations reflect electroweak gauge invariance, which we discussed previously. 

A useful comparison is to note that the Drell-Yan rates for selectron and stau pairs, 
shown in Table IX, are in the few tens of fb range for many SPS points. One might have 
hoped to use just the WBF sneutrino pair signals, which would often give missing energy 
plus forward-tagged jet signals, but because of the ultra-low signal rates, this will again be 
swamped by WBF Z ^ vv boson production [42]. 
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TABLE VI: Cross sections [fb] for WBF sneutrino pair production at LHC, for all MSSM bench- 
mark SPS points. Cross sections are shown to two significant digits or rounded to the nearest 
attobarn, and those smaller than half an attobarn are shown as zero. 



SPS 


la 


lb 


2 


3 


4 


5 


6 


7 


8 


9 




0.052 
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0.022 
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0.022 


0.007 


0.004 




0.045 


0.006 





0.021 


0.001 


0.017 


0.008 


0.069 


0.023 


0.001 


f+ff 


0.053 


0.016 





0.023 


0.006 


0.019 


0.008 


0.075 


0.025 


0.002 




0.007 


0.003 





0.001 


0.001 


0.002 





0.002 





0.001 


"^2^2 


0.043 


0.006 





0.014 


0.003 
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0.021 
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TABLE VII: Cross sections [fb] for WBF charged slepton pair production at LHC, for all MSSM 
benchmark SPS points. Cross sections are shown to two significant digits or rounded to the nearest 
attobarn, and those smaller than half an attobarn are shown as zero. 
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0.026 


0.003 





0.006 


0.001 


0.010 


0.002 


0.010 


0.003 


0.002 


fiVr 


0.005 


0.003 





0.001 


0.001 


0.001 





0.002 





0.001 




0.023 


0.003 





0.006 


0.001 


0.009 


0.002 


0.010 


0.003 


0.001 



TABLE VIII: Cross sections [fb] for WBF charged slepton plus sneutrino pair production at LHC, 
for all MSSM benchmark SPS points. Cross sections are shown to two significant digits or rounded 
to the nearest attobarn, and those smaller than half an attobarn are shown as zero. We neglect 
showing the charge-conjugate processes as the cross sections are trivially small. 
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SPS la 


SPSS 


process 


DY 


WBF 


DY 


WBF 




22.5 
29.0 


0.052 
0.045 


2.49 
14.3 


0.007 
0.023 




34.4 
18.3 


0.053 
0.043 


16.0 
2.40 


0.025 
0.006 



TABLE IX: Cross sections [fb] for Drell-Yan (DY) v. WBF slepton pair production at LHC, for two 
representative MSSM benchmark SPS points. Cross sections are shown to two significant digits or 
rounded to the nearest attobarn, and those smaUer than half an attobarn are shown as zero. The 
DY results were calculated with Prospino [5] at NLO. 

IV. CONCLUSIONS 

We have presented a thorough investigation of weak boson fusion production of colorless 
SUSY particles at the LHC. We find the cross sections to be almost universally unobservably 
small for all SPS scenarios, usually at the few-attobarn level, especially in the case of sleptons 
and sneutrinos. The smallness of the cross sections is partly due to large cancellations which 
take place at the amphtude level between WBF-type diagrams and bremsstrahlung diagrams. 

There remain two or three exceptions, where the rates are potentially interesting: X2X2 
production, which would give a highly distinctive four-lepton final state in many MSSM 
scenarios; and same-sign chargino production xtxt- The latter case is especially intriguing, 
because it can constitute definitive proof already at LHC that the ncutralinos in the t-channel 
are Majorana fermions. This is a crucial test for a candidate SUSY discovery. It might also 
provide information on the relative hierarchy of the wino and higgsino mass parameters. We 
sketched how such an analysis would proceed, identifying the dominant final state channels 
as a function of MSSM parameterization and outlining the major relevant backgrounds. 

To perform these calculations, we developed the matrix element generator SUSY- 
MADGRAPH which includes the complete MSSM with i?-parity and without additional CP 
violation, and does not assume any particular SUSY breaking scheme. It relies on the input 
from any SUSY spectrum generator in the SLHA format. We tested the MSSM implementa- 
tion by comparing with the literature for all known 2^2 collider scattering processes, and 
additionally via unitarity of scattering amplitudes at high energy and U{1)em gauge invari- 
ance. We furthermore derived analytical sum rules for neutralino and chargino electroweak 
gauge couplings from unitarity in all modes of VV — > xx scattering. 

We identified an issue of electroweak gauge invariance using SLHA input during develop- 
ment, on account of the unitarity checks: the default output of SUSY spectrum generators, 
combined with the default use of electroweak parameters for collider calculations, forms an 
inconsistent set of electroweak parameters. This misalignment can lead to sizable deviations 
in physical cross sections, relative to the known level of QCD perturbative uncertainty in 
the overall rates. SUSY-madgraph uses electroweak parameters taken from the neutralino 
and chargino mixing matrices to form a consistent set. 
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APPENDIX A: MADGRAPH II TECHNICAL DETAILS 

MADGRAPH [43] is a package which generates matrix elements for a user-input scattering 
process at a specified order in and a. The Fortran output it provides calls the HELAS [44] 
library of helicity amplitude subroutines, which can calculate any possible dimension-4 La- 
grangian term. It is extremely versatile, allowing as many external particles as one has 
computing power to handle, and allowing the user to require certain intermediate states 
while restricting others. It maintains all helicity correlations throughout. 

The original version had the Standard Model explicitly incorporated into the Fortran 
code. To facihtate additions beyond the Standard Model such as SUSY, MADGRAPH was 
rewritten such that it could read model information from an input file, and also handle 
Majorana fermions. These new features are incorporated into the release of MADGRAPH H , 
designed to work with MADEVENT [21], and briefiy described below. The new version 
also now provides sub- amplitudes in leading-iVc color fiows that allow it to be interfaced 
to standard parton-shower Monte Carlo programs such as pythia [45], HERWIG [46] or 

SHERPA [33]. 

The existence of Majorana fermions in SUSY required another significant modification 
to the code. MADGRAPH 's use of helas requires a well-defined continuous fermion fiow 
for the calculation of amplitudes and their proper interference. This continuous flow is not 
automatically satisfied by processes that involve Majorana fermions. Therefore, for a given 
process MADGRAPH H defines a continuous fermion fiow [47], chosen randomly from the 
two possible direction for any fermion line. For every fermion, MADGRAPH also defines a 
charge conjugate fermion with the opposite flow. Using this complete set of fermions, and 
requiring continuity of the flctitious fermion flow, MADGRAPH n is able to generate all of 
the appropriate diagrams, with the proper interference structure. In this scheme, conflicting 
fermion flows ("clashing arrows") are avoided by choosing a fermion flow randomly from 
the two possible directions, charge-conjugating the external wavefunction on one end of the 
fermion flow to match the chosen direction, and calculating the charge-conjugated vertex at 
the point of clashing arrows. 

The Standard Model Lagrangian in MADGRAPH H is specifled in two files, defining the 
particle content and the dimension-4 interactions, respectively. The user must set the defined 
couplings in a separate Fortran routine, with values to be passed to the matrix elements by 
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common block. The first file, particles.dat, defines all particles in the model. The next few 
lines show the format of the space-delineated file. 



d 


d~ 


F 


S 


ZERO 


ZERO 


T 


d 


1 


t 


t~ 


F 


S 


TMASS 


TWIDTH 


T 


t 


6 


e- 


e+ 


F 


S 


ZERO 


ZERO 


S 


e 


11 


g 


g 


V 


C 


ZERO 


ZERO 
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w- 


w+ 


V 


W 


WMASS 


WWIDTH 


s 


W 


-24 


h 


h 


S 


D 


HMASS 


HWIDTH 


s 


h 


25 



p uu~dd~ss~cc~g 



The first two entries define the characters to be used for the particle and antiparticle 
respectively. Notice in the case of a gluon it is the same symbol. The third entry gives 
the spin of the particle, choices are F for spin 1/2 fermions, V for spin 1 vectors and S 
for spin scalars. The fourth entry defines the type of line in the Feynman diagram, 
SCWD correspond to solid, curly, wavy, and dashed respectively. The next two entries are 
text strings representing the particles mass and width. The seventh entry gives the color 
information, valid options arc STO, representing color singlet, triplet or octet. The eighth 
entry is the text to be displayed on the Feynman diagram labeling the line, and the last 
entry is a number representing the particles ID to be used in interfacing with the QCD 
Les Houches Accord. In this format the complete set of Standard Model particles can be 
specified with just 17 lines. 

At the LHC it is particularly useful to generate a set of subprocesses by summing over 
groups of particles, e.g. partons in the proton, or final-state jets. The last line above defines 
the symbol p to represent the sum over the light quarks and giuon. There can not be any 
spaces between the particles that are to be summed over. 

The second file necessary to define a model is interactions.dat. As the name implies this 
file contains information about the allowed interactions, madgraph ii currently allows for 
3- and 4-particle vertices to be implemented using the space-delimited format shown below. 

d d g GG QCD 

g g g G QCD 

g g g g G G QCD QCD 

d d a GAD QED 

d u w- GWF QED 

u d w+ GWF QED 

The first line shows the vertex for a ddg QCD interaction. The first three entries are the 
particles appearing at the dimension-4 vertex. The ordering convention for fermion-fermion- 
boson interactions is first the incoming fermion, followed by the outgoing fermion, and the 
outgoing boson last. This is also illustrated in the interactions of the and W~ . The 
fourth entry contains the coupling strength for the interaction as it appears in the Fortran 
file, and the last entry specifies the type of coupling. This final string can be used to limit 
the diagrams to a certain order if a or a^. The four-gluon vertex illustrates a four-particle 
interaction. The first four entries define the particles involved in the interaction, the product 
of the fifth and sixth entries gives the couplings for the interaction, and the last two entries 
specify the couplings involved. 
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This is all the information required to define an interaction. The color structure of the 
vertex is inferred based on the color of the particles, and the Lorentz structure is inferred 
based on their spins. This format is compact, and the complete Standard Model (using 
diagonal CKM and MNS matrices) is specified with only 58 fines. 

SUSY-MADGRAPH is the MADGRAPH II package with two data files (as described above) 
which completely specify the MSSM Lagrangian. Moreover, it includes a Fortran routine 
which reads in the MSSM parameters from an SLHA files, and another which calculates the 
MSSM couplings. The input file may be generated by any MSSM spectrum generator. If 
one wants to include particle decays in the matrix elements, their widths must be specified 
in the SLHA file, and are read in automatically by the code. The couplings are assumed to 
conserve CP, but are written in the include file in such a way as to facilitate later inclusion 
of CP-violation, if the user chooses to add it. Chargino and neutralino mixing matrices 
are defined to be real. Negative neutralino masses appearing in the propagators do not 
introduce any errors for the Majorana scheme we implemented. 

APPENDIX B: NEUTRALINO AND CHARGINO FEYNMAN RULES 

In this appendix we give the explicit Feynman rules and couplings needed to evaluate 
the example sum rule for the process W~Z — > Xi)('i given as an example at the end of 
Section II B. Unfortunately, there are (at least) two sets of conventions for the chargino and 
neutralino mixing matrices present in the literature, namely Refs. [22, 23] and Ref. [24]. We 
give all formulas in both conventions. 

1. Neutralino and Chargino Mixing 

The chargino mass matrix according to Ref. [24] reads 

/ M2 V2mwC/3 \ 
Mc = I . (Bl) 

\ V2mwS/3 jJi J 

and can be diagonalized using two unitary matrices U^'^ MqU^ = diag(m~-). 

The unitary mixing matrices f/f^/j can be expressed in terms of two mixing angles. As 
long as these mixing matrices are real, one of the mass eigenvalues can in principle be 
negative. In that case a phase can be introduced in one of the mixing matrices {e.g. U^), or 
(as long as CP is conserved) we can perform the matrix element calculation with negative 
mass eigenvalues, simply by analytically continuing the expressions for the matrix element. 
Note that in HELASwe have to provide positive mass values for the spinor calculations. The 
analytic continuation holds only for the matrix element. 

The chargino mixing matrices can easily be translated into the conventions used in 
Refs. [22, 23], with the chargino mixing matrix Eq. (13). In that case we use the uni- 
tary matrices U and V to diagonalize the transpose of the chargino mixing matrix defined 
in Eq. (Bl) U* A4qV~^ . The translation rule becomes: 

C/f = C/t , Ui^V"" . (B2) 
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The neutralino mass matrix in both sets of conventions is identical, so we just copy the 
expression from Eq. (13): 



M 



N 



Ms 





\ mzSyjS/3 



mzCwCis 
-mzCwSfj 



-mzSyjCp 
mzCyjCp 




mzSwSp \ 
-mzCyjSp 





(B3) 



It is diagonahzed as U^^ M^U^ — diag(m^o) . Because the neutrahnos are Majorana 
fermions, the mass matrix Mjv is symmetric. Hence, up to a matrix of phase factors P, 
two unitary matrices and can be chosen as = Uj\fP* and = U^P. Again, if 
CP is conserved and we are wilhng to work with negative neutrahno mass eigenvalues using 
analytic continuation of the matrix element expression, we can choose P = 1. Alternatively, 
we can absorb the phase of the mass eigenvalue into a then-complex mixing matrix. 

Again, we can translate the neutralino mixing matrices in the conventions of Refs. [22, 23] . 
There, the mass matrix in the bino-wino basis is diagonahzed by a unitary transformation 
N*Aij^N~^ . The mixing matrices are now related by: 



N 



(B4) 



2. Couplings for interactions with gauge bosons 



Using the above definitions, we write out the couplings of gauge bosons to neutralinos and 
charginos for both set of conventions defined above. The symbol U without any superscript 
refers to the chargino mixing matrix according to Refs. [22, 23]. Note that the fermion order 
in HELASand in the sum rules is while in the SUSY-madgraph file interactions.dat 

it is 5'^"/°"*^; where /in,out are not necessarily both particle (as opposed to antiparticle). 
Instead, they are the particle or antiparticle for the rule as read directly off a Feynman 
diagram of the interaction vertex. This is important to remember in constructing chargino 
rules for SUSY-madgraph . The mixed chargino-neutralino-Vl^ couplings are: 



9l 



9r 



xjx°w 
9l 



x~x°w 
9r 



Ni2 Uji + -^Ni3 Uj2 

v2 



2j 



(B5) 



(B6) 



All neutralino and chargino mixing matrices are defined above and g is the usual weak gauge 
coupling. The chargino couplings to a photon or a Z boson are proportional to the gauge 
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couplings gz — g/cy,, with the usual weak mixing angle Cy, — cosOw and = sin^^^: 



Xi Xj z 

9l 



x-x-z 



9r 



„Xi Xi A 

yL 



9z 
9z 
9z 
9z 



Finally, the neutralino coupling to the Z is: 



xlx'z 1 



9l 



- -I9z 



k);k)3,-k);k) 



'4j 



-gz [iv,3 iv;3 - iv,4 iv;J , 



9r 



xlxV 



= -9l 



X^iX^.Z 



(B7) 



(B8) 



(B9) 



(BIO) 



(Bll) 



Due to the Majorana nature of the neutralinos, the indices i,j should be taken as i < j, 
and there is an additional factor of 2 for i = j in Eqs. (BIO) and (Bll). For the sum rules 
in Feynman gauge we need the Goldstone couplings to W,Z , neutralinos and charginos: 



9 



W-Zw+ 



-9 



W+Zw- 



-igzmws^ 



(B12) 



9l 



XiXj t^"* 



Xj x°^^ 



9r 



(B13) 



'2j 



C/3 



'7=2' 



n:,uj2 - V2N*,Uj, + '-^n:,uj2 



C/3 



XiXj t^"* 

9r 



= \9l 



Xj Xit^+ 



(B14) 



i 



(US),, (fg),, + V2 (i/f ); (US),. + (c/?), (US) 



'2j 



«/3 



S/3 
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